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ABSTRACT 

Consider  a finite  waiting  room  shared  by  several  servers.  Such 
a system  may  approximate,  for  instance,  a packet  switch  in  a 
commun i cati ons network  or  a buffer  shared  by  producer-consumer 
parallel  processes. 

It  is  assumed  that  a customer  is  lost  if  he  cannot  enter  the 
waiting  room.  Furthermore,  the  waiting  room  is  partitioned  into 
(L+l)  areas,  where  L is  the  number  of  servers:  each  of  the 
first  L areas  is  reserved  for  the  exclusive  use  of  one  server, 
the  (L  + l ) s t is  common  and  may  be  used  by  all  servers.  We  refer 
to  this  allocation  policy  as  ( L+l ) -s har i ng . 

The  steady-state  probabilities  are  determined  and  the  performance 
of  ( L + l ) -shari ng , in  terms  of  the  probability  that  a customer  is 
lost,  is  compared  to  three  other  policies: 

-restricted  sharing  (an  upper  bound  is  imposed  for  the  queue 
length  of  each  server), 

-unrestricted  sharing  (the  (L  + l ) s t area  occupies  the  whole 
waiting  room), 

-no  sharing  (the  (L  + l ) s t area  is  empty). 

It  appears  very  clearly  that  the  last  two  policies  are  not  the 
best;  a choice  between  ( L+l ) -shari ng  and  restricted  sharing 
would  depend  on  the  objectives  of  the  designer. 

Key  words: 

Exponential  queue,  finite  storage,  space  allocation  policy, 
numerical  comparisons 
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In  this  piper,  we  consider  the  problem  of  managing  a finite 
storage  space  used  concurrently  by  several  independent  servers. 
Systems  which  present  this  characteristic  appear  in  computer 
related  areas. 

-Consider  a packet-switched  communications  network  such  as  the 
ARPA  [8]  or  CIGALE  [10]  networks.  Each  packet  switch  has  a 
store  of  a finite  number  of  buffers,  the  buffers  are  shared 
by  all  the  output  links. 

-Parallel  processes  linked  in  producer-consumer  pairs  communicate 
via  a buffer  where  information  generated  by  the  producers  may 
queue  until  being  used  by  the  consumers. 

We  shall  assume  that  the  storage  space  is  made  up  of  M buffers; 
there  are  L output  processes.  Messages  are  produced  from  outside 
the  system,  destined  to  a specific  output  process  and  can  fit 
exactly  in  one  buffer.  We  assume  furthermore  that  if  an 
incoming  message  cannot  enter  a free  buffer,  it  Is  lost. 

Arrivals  of  messages  to  the  jth  output  process  form  a Poisson 
process  with  parameter  X.;  the  time  for  the  jth  output  process 

J 

to  service  one  message  is  exponential  u . ; all  random  variables 

J 

are  independent.  We  analyze,  the  system  i n • steady-sta>te  .- 
We  shall  defer  discussion  of  our  assumptions  and  of  their  effect 
on  our  conclusions  till  Section  6. 

One  has  to  define  a policy  to  allocate  the  storage  space  to 
the  different  output  processes.  We  shall  describe  first  the 
simplest  two  policies. 

Unres tri cted  sharing:  a message  may  join  the  system  whenever 


w 
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the  storage  space  is  not  filled;  no  other  restriction  is 
i mposed . 

No  sharing:  the  storage  space  is  partitioned  into  L disjoint 
areas;  a message  may  enter  the  system  only  if  the  area 
associated  to  its  output  process  is  not  filled. 

It  is  well  known  that  these  policies  may  lead  to  undesirable 
behavior  for  the  system.  Under  unrestricted  sharing,  one  of 
the  output  processes  may  monopolize  most  of  the  storage  space, 
if  it  is  very  slow  or  highly  utilized.  On  the  other  hand, 
under  no  sharing,  the  buffers  allocated  to  an  (almost)  inactive 
output  process  are  wasted:  they  are  not  used  by  their  process  and 
cannot  be  used  by  the  others. 

Two  policies  for  storage  allocation  have  been  defined  in  order 
to  reduce  the  impact  of  such  circumstances. 

Restricted  sharing  has  been  defined  by  Irland  in  1976  [6,7]: 
the  number  of  messages  in  the  storage  space  for  the  jth  output 
process  can  be  at  most  equal  to  K(<M).  Therefore,  if  k. 

J 

denotes  the  number  of  messages  in  queue  for  the  output  process  j, 
the  following  conditions  must  be  satisfied: 

L 

- z kr  . < -M-  

j = l J 


and  kj  < K,  for  1 < j < L . 

Irland  [6]  determines  the  steady-state  probabilities,  proposes 
algorithms  to  evaluate  the  normalizing  constant  and  the 
probability  that  a message  is  lost  and  gives  numerical  results. 
The  second  policy  has  been  defined  by  Dijkstra  in  1972  [4]. 

We  shall  refer  to  it  as  the  (L+1  )-sharinq  pol icy:  the  storage 
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space  is  partitioned  in  (L+1)  areas,  each  of  the  first  L areas 
is  reserved  for  one  output  process,  the  (L  + 1 ) s t area  is  common  and 
can  be  used  by  all  the  processes.  A message  may  enter  the  system  if 

a.  the  area  allocated  to  its  output  process  is  not  filled,  or 

b.  that  area  is  filled  and  the  common  area  is  not  filled. 
Therefore,  if  NQ  denotes  the  size  of  the  common  area  and  Nj  the 
size  of  the  area  associated  to  the  jth  process,  the  following 
conditions  must  be  satisfied: 

L 

Z N.  « M, 
j = 0 J 

kj  < N j + Nq,  1 < j <L  , 

kj  + kj'  - Nj  + Nj'  + V for  a11  0< j , j 1 <L , 


Z k.<  I N.  +Nn. 
j=l  J j=l  J U 


The  ( L+1 ) -shari ng  policy  is  used  commonly  in  the  literature  on 
parallel  processes  (see  for  instance  Devillers  and  Louchard 
[3]),  but  to  the  best  of  our  knowledge,  it  has  not  been 
analyzed  yet  in  the  context  of  stochastic  processes. 

If  L = 2,  the  restricted  sharing  and  ( L+1 ) -shari ng  policies  belong 
to  the  same  class:  if  N^N,,  and  K=N1+Nq,  they  are  identical. 

If  L>2,  those  two  policies  are  quite  different,  if  one  does 
not  consider  the  trivial  cases,  =N2= . • • =NL=K=M/L  , NQ=0 
(no  sharing)  and  N-j  =N2= . . . = NL=0  , K=NQ=M  (unrestricted  sharing). 
One  may  roughly  rank  the  four  policies,  in  order  of  increasing 
degree  of  sharing,  as  follows:  no  sharing,  ( L + 1 ) - s hari ng  , 
restricted  sharing  and  unrestricted  sharing. 
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In  Sections  1 to  3,  we  determine  the  steady-state  probabilities, 
which  always  exist  since  the  storage  space  is  finite  and  also 
the  probability  that  a message  is  lost,  we  analyze  the  special 
case  L=2  in  greater  detail. 

We  then  restrict  our  attention  to  balanced  ( L+l ) -s hari ng  policies: 
policies  for  which  N^=N23...=N^  (=N)  and  we  define  an  optimal 
solution  as  a pair  ( N Q , N ) which  minimizes  the  probability  that 
a message  is  lost. 

In  Section  5 we  present  some  of  the  numerical  results  we  have 
obtained  in  order  to  compare  the  performances  of  the  different 
policies.  It  appears,  as  was  noted  already  by  Irland  [6],  that 
the  unrestricted  sharing  and  no  sharing  policies  should  be 
omitted  from  further  practical  consideration.  In  particular, 
the  unrestricted  sharing  policy  proves  to  be  extremely  unstable 
under  heavy  traffic  conditions. 

Irland  defines  the  square  root  rule  which  is  a suboptimal 
restricted  sharing  policy,  i.e.  K=M//L.  It  is  interesting  as 
it  yields  good  performances  and  does  not  depend  on  the  system 
parameters . We  shall  introduce  in  Section  4 a similar  sub- 
optimal  balanced  ( L+l ) -shari ng  policy:  N=M/(L+/U),  Nq=M/(1+/T). 
The  relations  between  those  four  policies  (two  optimal  within 
their  classes,  two  suboptimal)  are  too  complex  to  be  summarized 
in  this  introduction. 

Finally  in  Section  6,  we  present  some  conclusions  and  we  discuss 
the  assumptions  we  have  made  in  our  model,  the  limitations  of 
the  present  approach  and  the  robustness  of  our  conclusions. 
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Steady-state  Probabilities 

The  state  of  the  system  is  represented  by  a vector  k_= 

( k.|  , l<2  . . . • »k|_ ) , where  is  the  number  of  messages  for  the  jth 
output  process  (in  short:  the  number  of  j-messages). 

Let  S represent  the  set  of  all  admissible  states: 

S={kJ OskjSNj+Ng  for  all  lsjsL, 

0<kj  + kj  , sNj+Nj  ,+N0  for  all  lsj,j‘<L; 

0<k  . + k . , +k  . „<N  . + N . , +N  . „+Nn  for  all  different 

j , j 1 »j  " ,1  <j  , j 1 , j "<L  ; 

L L 

0<  Z k.  < Z N.  = M}. 
j=l  J j=0  J 

Furthermore,  we  consider  the  probabilities 


P[k_,t]  = P[at  time  t,  the  system  is  in  state  Ic  ] , 

P[k_]  = lim  P[k,t],  and  the  quantities 
t *►« 


Xj/uj,  for  1 < j <L  . 


Since  the  storage  space  M is  finite,  the  steady-state 
probabilities  P[J<]  always  exist.  It  is  possible,  although 
tedious,  to  write  the  system  of  equilibrium  equations  which, 
together  with  the  normalizing  equation,  uniquely  determine  the 
P[kJ‘s.  This  system,  however,  is  not  particularly  enlightening 
and  we  shall  not  use  It  to  determine  the  steady-state 
probabilities  but  we  rather  use  the  local  balance 
equations . 

The  local  balance  equations  (see  [1])  express  the  condition  that 
in  steady-state,  the  rate  of  entrance  into  a state  of  the 
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system  by  arrival  of  a message  of  any  given  type  is  equal  to 
the  rate  of  exit  from  the  same  state  by  processing  a message 
of  the  same  type.  Formally 

x.  P[k{j)]  = Wj  P[k] 

for  all  j,  for  all  keS  such  that  k/^eS,  where  kj J ^ is  defined 
by 

(k_^J  ^ ) i = k.  , for  all  i/j, 

• kr'- 

This  system  has  the  solution 

1 L k . 

P[k]  = jr  n p,  J , 

L j = l 3 

where  C is  the  normalizing  constant  defined  by 

L k 

( 1 ) C = z n p . J . 

keS  j=l  J 

Let  it  and  it.  respectively  denote  the  steady-state  probabilities 

J 

that  a message  may  not  enter  the  system  and  that  a message  with 
destination  j may  not  enter  the  system.  One  has 
L 1 L 

(2)  TT  = ( z A.)  z A. IT., 

j=l  3 j=l  J J 


a *'£  ‘ P[k], 

3 ke  s 

J 


where  Sj  = { k_e  S | k ^ > N ^ and  there  exists  D<={1  ,2  , . . . ,N> 


such  that  Z k.  = z N.+Nn}, 

1 e D 1 ieD  1 0 

in  other  words,  S!j  is  the  set  of  states  such  that  all  the  buffers 
allocated  to  the  jth  output  process  and  all  the  common  buffers 
are  filled. 
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The  throughput  T is  defined  as  the  expected  number  of  messages 

that  enter  the  system  per  unit  of  time  in  steady-state,  i.e. 

L 

T = Z X . ( 1 - it  . ) - 7 is  clearly  given  by 

j = l 3 3 

L 

(3)  T = ( E XJd-ff). 
j=l  3 


2.  Special  case:  2 Output  Processes 


•v 


There  is  in  general  no  closed  form  expression  for  the  normalizing 
constant  C,  nor  for  the  loss  probabilities  it  . . If  L = 2 however, 

J 

one  can  easily  show  that 

-1  N,+N  +1  N2+N0+1 

(4)  C={  (1-p-j  ) (l-p2)  (p-|-p2)  > { (Pi -P2)  O "Pi  -p2  ) 

N,+l  N2+l  Nn  N 

+ Pi  P2  (p2  °(1-p2)-p-|  (1-Pi))}» 


-i  -i  -i  N2+Nn+1  -i 

(5)  it i = C P1  (l-p^)  (1-P2  )-°l  (1-P-J  ) > 

-i  -l  N, +Nn+1 

(6)  tt2  = C p2  ( 1 - p ! ) ( 1 - p 1 )-P2  ( 1 - p 2 ) , 


and 


-1 


(7)  ir=  1 - ( A .j  +X  2 ) (y^y.?) 


1-P 


W1 


+ CCCx1+x2 )}"  {;!-,  f— 


1-P 


W1 


+ y. 


1-P- 


Theorem 


For  a fixed  value  of  N2,  if  Ni  is  considered  as  a real  number 
in  [0,M-N2],  then  i as  a function  of  ^ has  a unique  minimum 
in  [0,M-N2]. 

Moreover,  if  for  some  value  N2  = N20,  it  is  minimum  for  N-^M-^, 
the  same  is  true  for  all  N2>N2°. 
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Proof : The  purely  technical  proof  of  this  theorem  is  presented 

in  Appendix  A. 

A completely  analogous  theorem  holds  for  » as  a function  of 
for  fixed  since  the  expression  for  it  is  symmetrical  in  the 
indexes  1 and  2.  These  theorems  and  other  technical  properties 
make  it  possible  to  determine  numerically  the  values  (N^.N^.N^) 
which  minimize  it  without  having  to  compare  all  the  possible 
sol utions . 

We  shall  not  present  detailed  numerical  results  for  this  special 
case  but  we  mention  that  it  appears  that 

(8)  lim  (maxiNg,^  ) > 0.5  if  p-]?*p2,c>1?sp2_1  , 

M-*-® 

where  N.  denotes  the  optimal  value  for  N^.  Moreover, 

maxtNg.N^  ,N2>  - NQ  if  p^p2<1, 

N.|  if  p i p 2 > 1 , p-|<p2, 

N2  i f p-|P2>1  , p 2<p-|  . 

The  rate  of  convergence  depends  on  the  ratio  y-j  / u2 . If  that 
ratio  is  close  to  1,  the  inequality  (8)  is  observed  for  values 
of  M as  small  as  10.  In  other  words,  we  have  an  asymptotic 
result  which  shows  that  for  the  optimal  solution,  one  of  the 
three  areas  occupies  at  least  one  half  of  the  buffer  space. 

3.  General  Case:  The  Loss  Probabilities  it. 

J 

For  L>2,  no  closed  form  expression  has  been  determined  for  the 
normalizing  constant  C or  the  loss  probabilities  ir . . We 

J 

rather  present  in  this  section  the  recurrence  relations  which 
may  be  used  to  compute  those  quantities  efficiently. 

Let  us  denote  by  S [ K ; N i , . . . , N K ; N Q ] , the  set  of  admissible 
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states  for  a system  with  K output  processes,  N . buffers 

vJ 

reserved  to  the  jth  process  and  NQ  common  buffers,  and  by 

C[K;N-j  , . . . ,Nk;Nq]  , the  normalizing  constant  for  this  system. 

Furthermore,  S ' [K;N^ , . . . ,N^;Nq]  is  the  set  of  admissible 

states  such  that  some  messages  may  not  enter  (some  of  the 

output  processes  are  "blocked"),  S 1 ={keS[K; . . . ;NQ] | 

3Dc(l  , 2 , . . . ,K>  , Z k . = Z N.  + Nn}  . 

jeD  J jeD  J U 

A1  so,  H[K;  N1,...,N|<;  Nq]  denotes  a sum  similar  to  the 
expression  (1)  for  C,  with  S and  L replaced  respectively  by 
S'  and  K; 


S j [ K ; N i , . . . , ; NQ],  j<K,  the  set  of  admissible  states  such 
that  no  message  to  the  jth  output  process  may  enter; 

G.[K;  N-j , . . . » N K ; N 0 ] , j<K,  a sum  similar  to  the  expression  (1) 
for  C,  with  S and  l replaced  by  S"  and  K respectively. 

We  obviously  have  that 


TT 


j 


G i [ L ; N1,...,Nl;  Nq] 
C [ L ; N1  ,. . . ,Nl;  NqJ 


for  0<j<L. 


Lemma  1 : 

-1  W1 

C [ 1 ; N-j ; Nq]=  ( 1 - p -j  ) ( 1 -p U ). 

For  K> 1 , 

C [ K+ 1 ; N1  ,N2,. . . ,Nk,  W;  Nq] 

■(1-pk+1  )-1  O-pjJ+i ) c CK ; N1,N2,...,Nk;  Nq] 
N0  • 

+ p“  + 1 2 pJ  + 1C[K;  Mn  ,N2 Nk;  N n - i ] . 
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Proof : 


The  first  relation  is  obvious. 


To  prove  the  second  relation  formally,  we  observe  that 

S [ K+ 1 ; N7  ,N2,. . . ,Nk,  W;  Nq] 

H-l 

= C (Jc.i ) I k_eS[K;  N1,...,Nk;  Nq]}] 


u[  U { ( k ,W+i ) | keS[K;  N,,...,N  ; Nn-1]}] 
1 = 0 “ ~ 1 K 0 

and  therefore 


C [ K+l ; N1  ,N2>. . . ,Nk,  Us  Nq] 

W-l  . K k 

= e . z n p,j 

1-0  K ' k«S[Ks. . . ;N0]  j=l  J 

N°  W+i  K ki 

+ e P„  j z n P,J. 

1=0  * 1 k«SfK;. . . ;Nn-1 ] j-1  J 


Informally  stated,  if  kK+1<W,  the  other  output  processes  are 
allowed  to  use  all  the  NQ  common  buffers,  if  W<kK+1sW+NQ,  the 
other  processes  are  allowed  to  use  ( N q - k K+ 1 + W ) common  buffers. 


Lemma  2: 


H [ 1 ; N] ; Nq]  = P] 


VM0 


for  K> 1 , 


H [ K+l  ; N 1 NK,W;N0]=(l-pK+1)'1(l-pjJ+1)H[K,N1  . . ,NK;NQ] 


+ PK+1  fgPK+lH^K;Nl  **  * * ’NK;N0~1-^ 
W+N  n K , N. 

j..  u n / 1 . \ " i / 1 . J \ 


+pk+i  n n-pj)'  n-pj  )• 
0 


■' 


1 1 


Proof : 

The  first  relation  is  obvious.  One  proves  the  second  one 
informally  as  follows.  The  first  term  on  the  right-hand  side 
corresponds  to  the  (K+l)st  output  process  using  less  than  W 
buffers,  all  the  Ng  common  buffers  being  used,  thereby  one  of 
the  other  processes  being  blocked.  The  second  term  corresponds 
to  the  (K+l)st  output  process  using  (*f+i)  buffers,  0<1<Ng,  all 
the  (Ng-i)  remaining  common  buffers  being  used  and  at  least 
one  of  the  other  processes  being  blocked.  The  last  term 
corresponds  to  the  (K+l)st  process  using  all  the  Ng  common 
buffers,  none  of  the  other  processes  being  blocked. 

Lemma  3: 

VNn 

g C i ; V n0]-Pi  1 u , 

for  j z2 , 

SjCJ ; N1,...,Nj_],  W;  Ng] 


W+N, 


C[j-1  ; N i , . . . ,N ; 0] 


Mq-1  . 

+ p i T P.-  H [ j - 1 ; N , , 
J i = 0 J 

for  K*j  , 

GjCK+1;  N1,...,Nk,  W;  Ng] 


,Nj-i ; v] • 


» Nn] 


+ PK+1  .^g  p K+l  Gj^K;  N1’--*’NK;  V1! 


Proof : 

The  first  relation  is  obvious.  The  second  can  be  informally 
proved  as  follows:  el ther  the  jth  process  uses  all  of  its 
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allocated  buffers  and  all  the  common  buffers,  the  other 
processes  may  use  any  number  of  their  allocated  buffers,  o_r 
the  jth  process  uses  all  of  its  allocated  buffers  but  not  all 
of  the  common  buffers  and  the  other  processes  use  the 
remaining  common  buffers. 

The  third  relation  is  proved  as  follows:  either  the  (K+l)st 
process  does  not  use  all  of  its  allocated  buffers,  the  other 
processes  use  NQ  common  buffers  and  the  jth  process  is 
blocked  or  the  (K+l)st  process  uses  at  least  all  of  its 
allocated  buffers,  the  other  processes  use  the  remaining 
common  buffers  and  the  jth  process  is  blocked. 

Remark : To  determine  C[K+1 ; , . . ;Ng]  , using  Lemma  1,  one  needs 
to  determine  C[K;...;0],  C[K;...;1],  ....  C[K;...;Ng]  and 
similarly  for  H and  G. 

Using  those  three  Lemmas,  it  is  possible  to  design  an  algorithm 
for  the  determination  of  the  probabilities  ir j and  therefore 
tne  probability  it  and  the  throughput  T.  The  algorithm  is 
presented  in  Appendix  B.  We  did  not  try  to  identify  the 
largest  problem  which  the  resulting  program  can  handle,  but 
a system  with  M=100,  L = 20,  NQ=20,  N-|  =N2=  • • . = N |_  = 4 was  solved  in 
about  0.5  seconds  of  CPU  time. 


Consider  the  problem  of  determining  a set  of  parameters 
L 

( N n , N , . . . ,N,  ) ( Z N.  = M)  which  maximize  the  throughput  T or, 
U 1 L j-0  J 

equivalently,  minimize  the  loss  probability  it.  Without 
further  analysis  of  it,  this  can  be  done  only  by  evaluating  tt 


1 3 

M + 1 

for  each  of  the  ( M ) possible  values  for  ( NQ  ,N^  , . . . ,NL ) . 

Such  small  systems  as  M=20,  L=5  begin  to  tax  the  computer,  even 
if  the  algorithm  avoids  repetition  of  lengthy  computations. 
Therefore,  we  restrict  our  attention  to  suboptimal  policies. 

A.  Optimal  balanced  ( L + 1 )- shari nq  policy 

An  optimal  balanced  policy  is  represented  by  (NQ,N)  and  is  the 
set  of  parameters  (Ng  ,N ^ , . . . ,N^  ) which  minimize  the  loss 
probability  tt  subject  to  the  following  constraints: 

n1=n2= . . .»nl-n, 

Nq+LN=M. 

The  determination  of  (NQ,N)  requires  the  evaluation  of  u for 
[M/LJ+1  systems. 

B.  Square-root  ( L+1 ) -shari nq  policy 

The  optimal  balanced  policy  depends  on  the  system  parameters  and 
requires  readjustment  as  traffic  character i sti cs  change.  As 
this  may  prove  difficult  to  do  in  practice,  one  would  set 
(Ng,N)  to  a fixed  value  which  yields  a reasonable  loss  probability 
it  over  a large  range  of  parameter  values.  We  find  such  a pair 
(Ng,N)  by  using  an  argument  similar  to  that  given  by  Irland  [6]. 
First  we  observe  that  losses  are  usually  insignificant  for 
small  values  of  the  p^.'s.  Therefore,  we  shall  analyze  the 
optimal  balanced  policy  in  the  special  case  where  p1=p2=.  . .=pl  = 1 . 
For  L=2,  one  can  show  from  Equations  (2)  and  (4)  to  (6),  by 
using  1 'Hospital ’s  rule,  that 

2 ( M-N  + 1 ) 

* " ( M+1 ) ( M+2 ) - 2N ( N+ 1 ) * 


V • • • 
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If  N can  take  any  real  value  in  [O',  ML"  ],  it  is  minimum  for 

N = M+l - /( M+ 1 ) (M+2)/2  , 
which  is  approximately  equal  to 

N - -S-  , 

2 + /2 


so  that  Nn=M-2N3 — — . 

0 1 + /2 

For  general  L>2,  we  define  the  square-root  ( L + l ) -s hari ng  policy 
as  follows: 


M or 

M 

U+/L  J 

, L+/L 

whichever  is  best  in  the  case  p . = 1 , l<j<L; 

J 

Ng=M- LN>0 . 

In  the  next  section,  we  examine  the  overall  quality  of  this 
policy.  It  seems  to  be  close  to  the  optimal  policy  if  p.  = l, 

J 

l<j<L,  although  we  do  not  have  a formal  proof  of  this  assertion. 
Let  us  mention  that  we  have  numerically  determined  the  optimal 
solution  for  the  following  systems: 

L = 3,4, ...  ,20; 

M = 20,30,. . . ,100; 

Aj  3 Wj,  1 < j <L . 

In  three  fourth  of  the  cases,  the  square-root  policy  is  optimal 
and  in  all  the  other  cases,  the  optimal  value  for  N is 


and  the  value  of  the  loss  probability  for  the  square-root  policy 
is  very  close  to  the  minimum. 
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Remark : Irland  [6]  determines  a square-root  restricted  sharing 

policy  which  is  close  to  the  optimal  for  p ^ = 1 , 1 < j < L : 

(9)  K = — . 

/L 

For  the  square-root  ( L+l ) -s har i ng  policy,  the  maximum  number  of 
buffers  that  any  output  process  may  use  is 

- M MM 

N + Nn  s — — + — — = — 

0 L + /U  1+/L  /L 

and  is  equal  (see  (9))  to  the  maximum  number  of  buffers  under 
the  square-root  restricted  sharing  policy.  Of  course,  under 
restricted  sharing,  two  or  more  output  processes  may  use 
simultaneously  M//T  buffers  each,  while  this  is  not  true  for 
( L+l ) -shari ng . 

5.  Numerical  Comparisons 

We  now  present  some  numerical  results  (see  also  Figures  1 to  14). 
Our  main  objective  is  to  compare  the  different  policies  we 
have  defined.  This  leaves  several  interesting  questions 
unanswered,  such  as  those  related  to  the  asymptotic  behavior 
for  M and  L-*®. 

As  six  policies  are  involved,  the  figures  tend  to  get  confused. 

To  alleviate  this  problem  somewhat,  we  adopt  the  following 
conventions. 

Conventions  : 

i.  The  loss  probabilities  are  presented  on  a logarithmic  scale. 

ii.  Results  for  the  unrestricted  sharing  policy  are  represented 
by  a "0", 


• • 
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-results  for  the  no  sharing  policy  are  represented  by  a 

-results  for  all  the  other  policies  are  represented  by  a 

continuous  line,  marked  by  a. 

"1"  for  the  optimal  balanced  ( L + l ) -shari ng  policy, 

"2"  for  the  square  root  ( L+l ) -s ha r i ng  policy, 

"3“  for  the  optimal  restricted  sharing  policy, 

"4"  for  the  square  root  restricted  sharing  policy. 

Moreover,  we  assume  in  all  of  the  examples  that  the  output 

processes  have  the  same  speed,  chosen  as  the  unit:  y.  = l,  1 < j < L . 

Hence,  p.aX.,  l<j<L. 
j j 

A.  Balanced  load  (figures  1 and  2) 

Consider  first  the  case  where  all  the  input  rates  are  equal: 
*1**2=’ ' •s*l*  Since  the  number  of  accessible  buffers  is  the 
same  for  each  output  link,  the  loss  probabilities  are  equal: 
ir  = n . , 1 < j < L . Figures  1 and  2 present  the  throughput  and  the 

J 

t loss  probability  for  L*3,  M= 1 0 and  for  different  values  for 

the  common  traffic  coefficient  p,. 

1 

For  small  values  of  (p-j<1),  the  throughput  for  the  no 
sharing  policy  is  substantially  smaller  than  for  the  other 
policies.  This  property  seems  rather  general  since  we  have 
systematically  observed  it  (see  figures  3,  7,  11)  whenever 
Pj<1 f IsjsL. 

For  large  values  of  p1  ( p ^ > 1 ) , the  unrestricted  sharing  policy 
offers  the  smallest  throughput,  followed  by  the  square  root 
restricted  sharing  policy.  This  suggests  that  the  unrestricted 
sharing  policy  is  the  most  sensitive  to  an  increase  of  congestion 
affecting  uniformly  all  the  output  processes;  and  that  the 
square  root  restricted  sharing  is  more  sensitive  to  such  a 

. • • a f 4 ^ - • • • • 4 - % - • • % *♦4  V 
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phenomenon  than  the  other  controlled  sharing  policies. 

B.  Unbalanced  load  (figures  3 to  10) 

For  the  next  example,  L=3  and  M=10  (figures  3 to  6)  and  20 
(figures  7 to  10).  Only  varies,  the  other  traffic 
coefficients  remain  constant  (p2=.5,  p3= . 6 ) . 

Notice  first  that  for  p^>l,  the  throughput  for  the  unrestricted 
sharing  policy  decreases  as  p1  increases.  This  property  seems 
to  be  a characteristic  of  the  situations  where  the  p . 1 s are  not 

J 

all  equal,  a justification  has  been  proposed  by  Irland  [6]. 

For  p i < 1 , the  global  loss  probabilities  n for  the  optimal 
balanced  ( L+l ) -sharing  policy  and  the  optimal  restricted  sharing 
policy  are  nearly  equal  (see  figures  4 and  8);  a similar  property 
holds  for  the  corresponding  square-root  policies.  Differences 
appear  for  p ^ > 1 (see  figures  3 and  7);  the  restricted  sharing 
policies  have  a larger  throughput  than  the  ( L+l )- shari ng 
policies.  However,  those  differences  get  smaller  when  M 
increases . 

In  this  case,  the  loss  probabilities  ir . are  not  equal;  figures 

J 

5,  6,  9 and  10  present  and  it ^ - Observe  that  for  all  the 
controlled  sharing  policies,  ir„  varies  on  the  interval  0<p-j<l 
much  less  than  , especially  so  for  the  square  root  policies; 
again,  this  property  is  general.  It  indicates  that  if  the 
traffic  is  stable  for  some  output  links  but  not  for  others, 
both  the  restricted  sharing  and  the  (L  + l )-sharing  policies 
reduce  the  effect  of  such  traffic  variations  on  the  stable  output 
processes.  Moreover,  one  observes  that  for  those  policies. 


For  our  last  example,  we  consider  a situation  where  most  of  the 
output  processes  have  a high  traffic  coefficient:  L=4,  M=20, 


p 2= . 7 , p . 8 , p^=.9,  p-j  varies. 

It  appears  first  (see  figure  11)  that  for  p-|>l,  the  square  root 
restricted  sharing  policy  does  not  perform  as  well  as  in  the 
preceding  example. 

Most  interesting  is  the  observation  one  can  make  from  figures  12 
to  14  ( p -j  < 1 ) : 

-the  global  loss  probability,  it,  takes  very  similar  values  for 
all  the  controlled  sharing  policies, 

-for  small  values  of  p^,  tt^  is  noticeably  smaller  under  (L  + l)- 
sharing  than  under  restricted  sharing, 

-there  is  less  difference  among  those  policies  for  the  other  loss 
probabi 1 i ti es . 

In  other  words,  if  most  of  the  output  processes  are  heavily 
loaded  and  a few  are  lightly  loaded,  the  latter  will  receive 
more  favorable  treatment  under  (L  + l )-sharing  than  under  restricted 
sharing  with  no  large  increase  of  loss  probability  for  the  other 
processes  nor  the  system  as  a whole. 

Again,  as  noted  earlier,  those  differences  decrease  for  larger  M. 
6.  Conclusions  and  Comments 

The  no  sharing  and  unrestricted  sharing  policies  prove  to  be 
completely  inadequate,  the  former  because  it  performs  poorly 
when  all  p.'s  are  smaller  than  1 (or  slightly  larger),  which 

J 

seems  a frequently  prevailing  condition,  the  latter  because  it 
presents  no  robustness  against  an  increase  of  traffic  for  some 
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or  all  of  the  output  processes. 

If  the  objective  is  to  minimize  the  global  loss  probability  t, 
the  optimal  restricted  sharing  policy  proves  to  be  best  by  a 
small  margin  in  the  majority  of  cases:  we  have  found  very  few 
(usually  unrealistic)  situations  where  the  optimal  balanced 
( L+l ) -s hari ng  policy  is  better. 

If  the  objective  is  to  have  a small  global  loss  probability  *• 
and  a small  loss  probability  for  lightly  loaded  output 
processes,  then  there  is  no  uniformly  best  policy.  In  cases 
similar  to  exampl e B,  restricted  sharing  gives  best  results, 
for  situations  similar  to  example  C,  it  is  ( L + l ) -shari ng . 

The  square  root  policies  have  the  advantage,  as  we  mentioned 
earlier,  of  being  independent  of  the  parameters.  They  perform 
well  in  general  and  there  does  not  seem  to  be  any  general  answer  to 
the  question  of  which  of  the  two  square  root  policies  minimizes  the 
global  loss  probability  it.  However  the  square  root  ( L + l ) -s  hari  ng 
policy  may  prove  to  be  preferable  if  one  takes  into  consideration  the 
risk  of  coalitions  (defined  informally  below):  coalitions  are 
impossible  under  square  root  ( L+l ) -s hari ng  if  M>L+/T,  but  may 
occur  under  square  root  restricted  sharing  if  L>3.  A coalition 
occurs  when  the  output  processes  have  a behavior  which  allows  a 
subset  of  them  to  always  use  all  of  the  M buffers,  preventing 
the  other  processes  from  ever  getting  any  message.  This  is 
reflected  to  some  degree  by  the  values  of  for  small  values  of 
p.j  in  example  C (figure  13). 

We  now  comment  on  some  of  the  assumptions  we  have  made  in  defining 
our  model . 
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Firstly,  we  have  assumed  that  arrivals  form  a Poisson  process 

- 

and  services  are  exponential.  To  the  best  of  our  knowledge, 
only  Fisher  [5]  considers  a general  service  time  distribution, 

L=2.  Fisher  proposes  an  approximation  for  the  steady-state 
probability  distribution  under  unrestricted  sharing  and  makes 
comparisons  with  the  no  sharing  policy  in  the  exponential  case. 
Extending  the  analysis  of  our  model  to  general  service  time  and/or 
interarrival  time  distributions  seems  to  be  a major  task;  before 
engaging  in  it,  some  measurement  of  the  robustness  (or  lack  of 
robustness)  of  the  results  for  the  Markovian  case  might  be 
appropri ate. 

Secondly,  we  have  considered  the  system  in  stationary  state. 

» 

This  may  not  be  appropriate  since  the  conditions  are  likely  to 
change  over  time  for  the  applications  described  in  the  intro- 
duction. However,  the  results  we  have  obtained  suggest  that 
the  square  root  policies  are  generally  close  to  the  optimal  and 
should  be  good  ones  provided  that  abrupt  changes  do  not  occur 
too  frequently,  especially  if  one  takes  into  consideration  the 
costs  of  adjusting  the  allocation  policy  to  a changing 
envi ronment. 

Two  of  our  assumptions  have  a crucial  influence  on  our  analysis. 

We  have  assumed  that  a message  that  cannot  join  the  system  is 
lost.  This  may  be  legitimate  for  a packet-switch  communications 
network:  such  a message  would  either  be  resubmitted  later  (and 

appear  to  be  a new  one)  or  transmitted  through  another  path  in 
the  network  (and  be  effectively  lost  for  the  system  under 
j consideration).  For  producer-consumer  parallel  processes,  this 

• . • I • «♦'  - - • • • • ••  * ' • •*»- 
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assumption  is  at  best  an  approximation.  A producer  that  cannot 
find  a buffer  to  enter  a message  will  wait  until  it  gets  a free 
buffer.  This  should  affect  the  steady-state  probabilities  but 
not  necessarily  our  qualitative  comparisons  of  the  different 
policies. 

We  have  also  assumed  that  all  random  variables  are  independent. 
This  affects  our  qualitative  conclusions.  For  instance,  this 
assumption  reduces  the  risk  of  coalition.  It  is  possible  to 
define  arrival  processes  which  are  not  independent,  such  that 
the  square  root  restricted  sharing  policy  and  even  the  optimal 
restricted  sharing  policy  (based  on  the  Markovian  assumption) 
perform  disastrously:  coalitions  are  certain  to  occur  and  the 
output  process  with  the  smallest  traffic  coefficient  never  gets 
any  messages. 

Furthermore,  situations,  not  covered  by  our  numerical  analysis, 
are  worthy  of  future  investigation.  For  instance,  would  our 
qualitative  conclusions  hold  if  the  service  rates  y . 1 s were 

J 

not  all  equal?  More  importantly  perhaps,  one  should  try  to  get 
asymptotic  results  for  large  values  of  M and  L;  intuitively,  it 
seems  that  for  fixed  L and  increasing  M,  the  differences  between 
policies  would  reduce.  This  is  supported  by  our  numerical 
results . 

Finally,  one  may  consider  that  the  messages  have  not  all  the 
same  importance  and,  for  instance,  assign  costs  for  loss  of  a 
message,  which  are  different  according  to  the  destination.  The 
function  to  minimize  in  this  case  is  the  expected  cost  per  unit 
of  time: 
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where  Yj  is  the  cost  for  loss  of  a message  with  destination 

j • 

In  [9]  we  analyze  the  no  sharing  policy  in  a small  system: 

L = 2,  exponential  services.  For  general  L and  ( L + l ) -s hari ng 
policy,  f can  be  determined  as  easily  as  it  (Section  3).  If 
costs  i.re  different,  one  cannot  restrict  the  search  for  an 
optimum  among  the  balanced  policies;  we  have  already  indicated 
in  Section  4 that  an  exhaustive  comparison  of  all  the  possible 
values  for  (Ng,N^,...,N,)  is  prohibitively  costly  except  for 
very  small  systems.  Therefore,  r would  have  to  be  analyzed 
further  in  order  to  determine  an  efficient  algorithm. 


23 


Appendix  A 

1.  Let  L=2.  Consider  N2  fixed  and  let  take  any  real  value 
in  [0,  M-N2].  it  as  a function  of  has  a unique  minimum  in 

[0,  m-n2]. 

Proof : 

As  N2+N-|  + N0  = M,  the  expression  (7)  for  it  can  be  written  as 
ir=Y+6g -J  ( N]  ; N2.  M ) / g 2 ( N 1 ; N2,  M) 

where 

y=1-(a^+X2)"^  (tJ-]+U2^  1S  independent  of  , 

5=(A^+A2)-1  is  positive  and  independent  of  , 

M-N, +1  M-N2+l 

g i = u i ( i - p 1 ) ( l -p  2 )+u2  0 -p2  ) ( 1 -p-|  )> 

_1  M-N2  + l M-N-j  +1 

g2=(pl"p2^  f ( P 1 -p 2 ) O -p 1 -P2  ) 

N , + 1 M-N, +1  M-N-+1  N?+l 

+ p^  P2  (1-P2)-P1  P2  (1-q-j  )}- 

Replacing  N ^ by  x,  it  is  easy  to  check  that 
tt1  (x)  = (h1  (x)  )2h  (x) 

where  h -j  ( x ) is  a function  of  no  interest  to  us  and  h(x)  is 
given  by 

_1  M-N2+l 

( A . 1 ) h(x)  = (p-|-p2)  (-y -|  ( 1 -p-|  ) (p -| -p2  ) p-|  log  p2 

2 N2+l  M-N2+l 

-u1 ( 1-p1 ) p2  p1  log  p2 

m-n2+i 

-u2 ( 1 _P2 H p-j -p2 ) ( 1 -p-|  )log  p2 

+y  ^ ( 1 -p-|  ) ( 1 -p2  )p-j  X+1  p2M  x+^log  p2 

m v+i  M~N9+1 

-[y1d-p1)(i-p2  x )+v2n-p2)n-pi  )] 

(i-p2)p1x+1nogpriogp2)}. 
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The  derivative  of  h(x)  is 

X+  1 - ] 

h'(x)  = -p1  (logp1-logp2)  (p1-p2)  x 
{w-j  (1-p2)(1-p2M  x+^  ) ( 1 -p  i ) 1 ogp  i 

2 m-n2+i 

+ u2 ( 1 -p2 ) ( 1 -p1  ) 1 ogp1  } 

and  is  positive  since  (l-pa)logp<0,  a>0;  hence  h(x)  is 
increasing  and  it  Is  trivial  to  conclude  that 
ir  is  minimum  for  N^=0  if  h(0)>0, 

ir  has  a unique  minimum  in  (0,  M-N2)  if  h ( 0 ) <0<h  ( M - N 2 ) , 
ir  is  minimum  for  N -j  = M- N 2 if  h(M-N2)<0. 

2.  If  h(M-N2°)<0,  then  h(M-N2)<0  for  all  N2>N2°. 

Proof : 

From  (A.l),  h(M-N2)<0  if  and  only  if 

N?+l  , M-N?+l 

(A. 2)  u1  (1-p1  )(1-p2  )(p1-p2)“  p1  C ( 1 - p -j  )logp2-(l-p2)logp1  ] 

M - N 2 + "I  _1 

-u2(1-p2)(1-p1  )(p1-p2)’  [ (p-j -p2  ) logp2+ 

m-n2+i 

(1-p2)p1  (logPl-logp2)]<0. 

As  (l-p)’^logp  is  an  increasing  function  of  p,  it  is  easy  to 

verify  that  the  coefficient  of  y-j  in  (A. 2)  is  negative. 

The  coefficient  of  (-y2)  is  negative;  it  can  be  written  as 
, M-N,+l 

(1-p-j  )’  (1-P-,  ) g ( N 2 ) where 

g ( x ) = ( 1 - p i ) ( 1 - p 2 ) ( p i ~ p 2 ) [ ( p i ■ p 2 ) 1 09p 2 + 

(1-p2)p1M-x+1 ( 1 ogp i -1 ogp  2 ) ] . 
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As  g ' (x)=-( 1-p1 ) ( 1-p2)2 (pt -o2)"1p1M’X+1 logp1 (logp] -logp2)  is 

positive  and  g ( M ) = ( 1 -p ] ) 2 ( 1 -o2 ) 2 ( p 1 -p2 ) ' 1 [ ( 1 -p 1 ) -1 p 1 1 ogp 1 - 
(l-P2)-1P2logp2]  is  negative  since  (l-p)_1plogp  is  a decreasing 
function  of  p,  g(x)  is  negative  for  all  x<M. 

Hence,  h(M-N2)<0  if  and  only  if  y-|/u2sb(N2),  where 

b(x)  = { (1  - p-j  ) ( 1-p2x+1  )p-jM  X+1  [ ( 1 -P] ) 1 ogp2- ( 1 -p2  ) 1 ogpi  ]}  ^x 

{(1-P2)(1-PiM  X+1)[(P|-P2)1°9P2  + ^1"p2^p1M"X  + 1^09p1"1°9p2^^-® 
After  tedious  algebraic  manipulations,  one  can  write  b'(x)  as 
b,(x)  = p1M"x+1{(l-p1)4(l-p2)2(prp2)'1[(l-p1)‘1logpr 

(l-P2)'1logp2]}‘1x{.(l.p2)(l-p2x+1)(l-p1)logp1g(2x-M-l) 

-(l-P1)(l-P1M'X+1)P2X  + 1n-P2)1°9P29(x)} 


which  is  clearly  negative  for  all  xsM.  Therefore,  b(x)  is 
decreasing  and  b(N2°)<y-j/y2  implies  b ( N2  ) -j  / u 2 for  a11  N2>N2°, 
which  completes  the  proof. 
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Appendix  B 

Denote  by  x.  an  Infinite  sequence  ( Xg  ,x-|  ,x2 , . . . ) and  let  the 
operator  * denote  the  convolution: 

j 

(x*z)j  = 2 Xi  yj-1  0<j. 

J i = 0 J 1 

The  Lemmas  1 to  3 (Section  3)  can  be  written  as  follows,  using 
Nk+i  instead  of  W.  To  simplify  notations,  we  use  CK(Ng) 
instead  of  C[K;  N-j  ,N2 , . . . ,NK ; NQ]  and  we  let  represent  the 
sequence  (C^(0),  CK(1),  ...),  similarly  for  H and  G.  Moreover, 

P.  j " ( 1 » P j » ^ j * •••)> 

and  we  define  the  sequence  v.  as  follows: 

j 


-j 

Lemma  1 ' 


v , ( n ) = z 


i = 0 


1 < j < L , 0<n. 


C,  (Nq)  = v^N^Ng)  0<Nq, 

NK 

CK  = vk(Nk-1)  CK-1  + PK  (p^C^.,)  2<K<L. 

Lemma  2' 

Ni+Nn 

H-j  ( NQ ) = Pl  1 U 0<Ng, 

nk< 


-K  = -K(NK_1)  -K-l  + PK  K*— K- 1 ^ 

N.  K-l 

+ P|(  £k  n — j ( N j “ 1 ) • 2<K<L  . 

J 1 


Lemma  3' 


ii.i  = Jtv 


— K , K = PK  (— K- 1 (°)'-K-l  (°))fi-K 


+ P|(  K(£.K*tL|(..-|  ) > 2<K<L 


— j , K = -K(NK_1)  — j , K-  1 + p K K^K*-j  ,K-1  ^ * 1-j<K-L 


To  compute  the  first  N components  of  £=x*£_  where  p_=(l,p,p' 
we  use  the  following  algorithm  (Buzen  [2]): 


CONVOL  (y,x,p,N); 


real  array  y ,x ; real  p ; i nteqer  N ; 


begin  i nteqer  k ; 


y(0) :=x(0) ; 


for  k : =1  step  1 unti 1 N do 
y(k) :=py(k-l )+x(k) ; 


This  algorithm  requires  N multiplications  and  N additions. 


Therefore,  to  determine  C and  G.  , , 1 < j < L , one  has  to  perform 

N 0 ( L - 1 ) ( L+7 )+M+2L  additions 
and  Nq(L-1 ) (3L/2+5)+M+L  multiplications. 

(Details  are  indicated  in  Table  1).  Moreover,  memory  require- 


ments are  not  very  large  if  one  discards  intermediate  results 


when  they  are  no  longer  needed. 
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C i (£-1  ar,d  ^-1 ) 


C„  (2sK<L) 


(£k  and  vK) 
convol uti on 
total 


H (2<K<L) 

convol uti on 
last  term 

tota  1 


Gk>k  ( 2 < K < L ) 


G . „ (2<K<L) 

convolution 

total 

total  all  j 


H1  and  Gl,l 


28 

addi ti ons 


VV1 


W1 


V3No+1 


3N, 


V1 


2N 


0 


2 ( K- 1 ) N 


0 


Table  1 


multipli cations 


VN0 


VN0 


V4N0 


4N, 


V 


3Nq 

3 ( K - 1 ) N, 


29 
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Figure  2 Loss  probability  it;  balanced  queue,  L = 3,  M=10 
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Figure  3 Throughput;  unbalanced  queue,  L=3,  M=10,  p2=.5,  p33.6 
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Figure  6 Loss  probability  for  2-messages:  unbalanced  queue; 

L=3,  M=10,  o 2 = • 5 , p 2=  • 6 
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Figure  10  Loss  probability  for  2-messages:  ir?;  unbalanced  queue; 
L=  3 , M=  2 0 , p 2= . 5 , p . 6 
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Figure  14  Loss  probability  for  3-messages:  tt  3 ; 

highly  unbalanced  queue;  J 

L = 4 , M=  20  , p 2 = • 7 i p 2= . 8 , p^=.9 
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